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164000-2-T 


Second Quarterly Report 

STUDY ON SPECTRAL/RADIOMETRIC CHARACTERISTICS 
OF THE THEMATIC MAPPER FOR LAND USE APPLICATIONS 


1. Objective 

The objective of this investigation is to quantify the performance 
of the TM as manifested by the quality of its image data in order to 
suggest improvements in data production and to assess the effects of 
the data quality on its utility forland resources applications. Three 
categories of this analysis are: a) radiometric effects, b) spatial 

effects, and c) geometric effects, with emphasis on radiometric effects. 


2. Tasks 

Four tasks have been established to address the above objective. 

The first three are to study radiometric performance, spatial perfor- 
mance, and geometric performance, respectively , while the fourth is to 
study spectral characteristics. In keeping with the identified objective, 
the radiometric performance study is the major task. 


3. Status and Technical Progress 

Radiometrically and geometrically corrected Thematic Mapper digital 
data (CCT-PT) and image products were received for two scenes during 
this quarter; Scene 40145-14492 (8 Dec Cape Cod), and Scene 40122-16234 
(15 Nov Oklahoma). The radiometrically , but not geometrically, cor- 
rected data (CCT-AT) have been ordered but not yet received. Raw data 
(CCT-BT) for these two scenes have not been released as of this time. 

William Mali la and Michael Metzler presented preliminary findings 
at the Second Landsat-4 Investigators 1 Workshop 11-12 January 1983. 

They also attended the Landsat-4 Scientific Characterization of Early 
Results Symposium, where the paper "Scan-Angle and Detector Effects in 
Thematic Mapper Radiometry" by Metzler and Mali la was presented, The 
paper describes the analysis and findings of this investigation to date; 
it is included in this report as Appendix A. 

A paper entitled "The Enhanced Earth Observation Potential Realized 
With the Thematic Mapper" was presented at the Landsat-4 Session of the 
1983 ACSM-ASP Annual Meeting on 17 March 1983, 
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3.1 Problems 

No significant problems were encountered this quarter. 


3.2 Accomplishments 

The accomplishments to date are described in detail in Appendix A. 
The findings of this reporting period are summarized here. 


3.3 Significant Results 

Scan-Direction and Scan-Angle Effects 

The observations made during the previous quarter of the change 
in mean signal level as a function of both scan-angle and scan-direction 
were studied further. The overall scan-angle effect was found to corre- 
spond to that expected based on atmospheric modeling and scene charac- 
terlstics. An initial, empirical correction model employing exponential 
decay was developed for the reflective bands in the two scenes analyzed 
(40049-16262 and 40037-16031). Band 6 was observed to have a significant 
scan-direction effect which was markedly different than that found in the 
reflective bands. 

Later-Detector Calibration 


A low frequency (period of several scan cycles) noise was discovered. 
This noise was most pronounced in Band 1, Detectors 4, 12, 10, and 8, 
having amplitudes of approximately 2,0, 1.5, 1.0, and 0.75 quantizing 
levels, respectively. This low frequency variation in mean signal ampli- 
tude was highly correlated among those four Band 1 detectors. Low fre- 
quency noise was also observed in Band 7, Detector 7; Band 2, Detector 1; 
Band 3, Detectors 1 and 16; and Band 5, Detector 10. For the non-Band 1 
detectors, the amplitude of the noise is £0.5 quantizing levels. This 
noise may be a result of interaction of the DC-restore process with 
coherent (high frequency) noise found in the system. 


3,4 Publications 
None 
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3.5 Recommendations 

It is recommended that investigations into the scan-direction 
effect ahd low-frequency noise be pursued. Any corrections resulting 
from such an investigation should be incorporated into a single radio- 
metric correction to minimize quantization and round-off effects. 


3.6 Funds Expended 


3.7 Data Utility 

Thematic Mapper has been shown to provide greater spatial and 
spectral resolution than the MSS. Correction of the small but signifi- 
cant effects described in this report will be helpful in many applica- 
tions of the data, essential in some. For studies of waterbodies in 
particular, where small signal changes can be significant and where 
Band 1 is very important, the low frequency noise and scan direction 
effects observed in Band 1 will reduce the utility of the data if they 
remain uncorrected. 
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SCAN ANGLE AND DETECTOR EFFECTS IN 
THEMATIC MAPPER RADIOMETRY 


Abstract 


The effects on Thematic Mapper (TM) radlometry of scan angle and 
Inter-detector differences were analyzed. Radiometric corrections 
currently performed were found to improve overall consistency of data, but 
some residual striping remains in the corrected data due to the 
quantization of signal values and other effects. A new type of banding 
was discovered which is related to the bidirectional scanning of TM. An 
initial empirical mode 1 / was developed for correcting this effect in Band 
1, but should receive additional development. The scan angle effects 
observed corresponded to those expected based on atmospheric 
considerations and scene characteristics. Low frequency scan-to-scan 
noise was detected in Band 1, Band 7 and to a lesser extent in Bands 2 and 
3. The Band 1 detectors which exhibited this noise showed strong 
correlation in their variation. 


INTRODUCTION 


This investigation was directed at quantifying and understanding the 
scan-angle and detector effects observed in TM radiometry. 

DATA AMD METHODS OF ANALYSIS 


The data set used for this study consisted of computer compatible 
tapes (CCT's) of raw data (CCT-BT), radiometrically corrected data 
(CCT-AT), and geometrically corrected data (CCT-PT) for two scenes, 
40049-16262 (North Central Iowa) and 40037-16031 (Arkansas). Major 
emphasis was given to the unresampled (CCT-BT and CCT-AT) data to allow 
analysis of individual scan and detector effects. Although analysis was 
concentrated on the Iowa scene because of its cloud-free, relatively 
homogeneous scene content, all analyses were performed for both scenes for 


confirmation of results. The analyses were performed over the full frame 
of data for each band, for each detector within each band and for each 
scan direction. 

Histograms of the frequency of occurrence of each quantized signal 
level were produced and analyzed for each step 1r. the correction process 
(CCT-BT to CCT-AT to CCT-PT). Signal means and variances were computed, 
again before and after calibration. Across track profiles of the mean 
signal as a function of pixel number were generated for each scan 
direction by dividing the scene Into forward and reverse scans and then 
computing the mean signal for each set of scans In a window 16 samples 
wide by 2992 lines (i.e., all forward or all reverse lines) long. 
Stepping the window across the scene produced relatively noise-free 
average scan lines exhibiting only gross scene content effects. Along 
track profiles were generated for each detector by computing the mean 
signal level for all the pixels in each scan line and then plotting the 
scan-line means for a given detector as a function of scan number (374 
scans x 16 detectors = 5984 lines In full scene, CCT-BT or CCT-AT). 


RESULTS 

* 

The overall quality of TM imagery was observed to be good. However, 
several effects were noted which could be significant in analysis of data 
from different scenes, from different areas within the same scene or for 
applications where small signal differences are important. 

Scan Angle and Scan Direction Effect s 

Scan angle effects caused by atmospheric backscattering, changing 
optical path lengths, and ground bidirectional reflectance effects were 
expected. Figure 1 illustrates the mean signal level as a function of 
scan angle (pixel number) for each of the seven TM bands. The sharp peaks 
present are primarily due to urban areas (at pixels 800, 3300), larger 
bodies of water (pixels 4200-5200), or numerous smaller bodies of water. 
The feature of interest, however, is the overall decline in mean signal 
level in going from the western side of the scene to the eastern side. 
The effect is greatest in Bands 1 and 2 which decrease approximately 10%. 
Band 3 exhibits a somewhat less apparent angular effect, and the remaining 
reflective bards demonstrate an even lesser effect and more influence of 
scene composition. The thermal band. Band 6, also exhibits a similar 
pattern, with larger responses at the western edge, presumably due to 
viewing a greater percentage of sunlit scene components. However, this 
full-frame average scan line for the thermal band masks a distinctive 
scan-direction effect which is discussed in detail later. 

Previous modeling work at ERIM for analyses of spaceborne scanner 
data have predicted atmospheric scan angle effects similar to those seen 
with this TM data [1,2]. Figure 2 is a plot of model results derived by 
us from the data set of Dave [3] for the MSS; their approximate TM 
equivalents are given. At the longer wavelengths, e.g., TM Bands 4, 5 and 
7, atmospheric scan angle effects are reduced, and canopy shadowing and 







bidirectional reflectance effects become the predominant drivers of the 
observed scan-angle effect. 

Quantification of the scan-angle effects in TM radlometry was 
complicated by a new effect which was discovered. Figure 3 illustrates 
the mean signal level as a function of scan angle for Bands 1 and 6 as In 
Figure 1, but In Figure 3 the forward and reverse scans are plotted 
'separately. While the dominant effect is the expected scan angle effect, 
a systematic droop of signal values during the active scan tends to 
Increase the observed scan angle effect for forward scans, and decrease it 
during reverse scans. This scan direction effect was detected to some 
extent for all of the reflective bands, and was most pronounced In Band 1. 
Band 6 also demonstrated a pronounced scan direction effect, but of a 
distinctly different nature (see Figure 3). 

In order to characterize the scan direction effect, the ratio of 
forward scan signal level to reverse scan signal level was computed as a 
function of scan angle (pixel). These data (Figure 4) were than fit to an 
exponential decay model of the form 



1 + A[exp(-KM) - 1] 


where 

p = Pixel number, counting from West edge of scene 

Pf»P r = Pixel offsets to convert "p" to relative minor frame number 

M = p>p for forward scans 

p p -p for reverse scans 

S Q (p) = Signal returned to sensor for pixel "p" 

S(p) = Corrected signal value for pixel "p" 

K = Time constant (reciprocal of the number of Dixols of active 

scan required for 63£ of the decay to occur) 

A = Factor determining magnitude of total decay 

The model was fit with the constraint that A and K be identical for 
forward and reverse scans. The result of fitting this model to Band 1 of 
the Iowa scene is illustrated in Figure 4, and in Figure 5 the predicted 
signal decay during forward and reverse scans is shown. From this model, 
one can see that the decay begins prior to the first pixel, having decayed 
approximately 1£ by that time. The total decay is approximately 2 . 25% by 
the end of the active scan. This results in the mean signal level of the 
forward scan being approximately 0.75 counts higher than the mean signal 
level of the reverse scan at the West end of the scene, and approximately 
0.75 counts lower at the East end of the scene for Band 1. 


For this case {Band 1, scene 40049-16262) the p- wnd p required for 
the fit resulted in the exponential decay beginning at tne time of DC 
restore (different times for the two scan directions). During DC restore, 
a dark (no radiance) level is presented to the detectors, and a DC level 
capacitor is charged to a value which results in approximately two signal 
counts. This charge is t<? be held for the entire scan, having a time 
constant on th.e order of minutes [4]. The, empirically derived 
scan-direction effect model suggests that this charge may be decaying much 
faster (time constant of the order of 10 msec). Two problems exist with 
this explanation: 1) 10 msec is much, much less than minutes and 2) the 
magnitude of the effect appears to be proportional to the mean signal 
level, while the DC restore capacitor simply adds a set signal level. 
Final explanation of the effect awaits further investigation. 

Any corrections for this scan-direction effect should be made an 
integral part of a single radiometric correction process to minimize 
quantization and roundoff errors. 

Inter-Detector Calibration Effects 


Horizontal striping is commonly found in imagery from sensors that 
scan arrays of detectors. To minimize the potential for striping among 
the 16 detectors in each band (four in Band 6), a histogram normalization 
technique is included in the radiometric correction processing. This 
technique attempts to equalize the means and variances of the quantization 
level histograms of each detector within a band by applying a gain and 
offset to the signals from each detector. Since this gain and offset is 
applied to a discrete function (hi stogram) , and the result is still 
integer, the number of occupied quantization levels is not changed by the 
correction; only the total range they span is changed. This requires that 
some quantization levels within the expanded range be empty. The effect 
of this correction is illustrated in Figure 6. Although the histogram for 
an individual detector acquires some empty bins in the correction process 
(right-center), the histogram of the entire band is noticeably smoother 
(left-center). Of course, the cubic convolution resampling done in the 
geometric correction processing tends to equalize the histogram bin sizes 
even more (left-bottom). 

The net effect of tNse corrections for a band with a fairly large 
variance, such as Band 4 , illustrated here, is virtual elimination of any 
detector striping in the radiometrically corrected data. On the other 
hand, if a band has a lower variance, the histogram normalization is 
unable to equalize the means and variances of the detectors because of 
quantization limitations, and low level residual striping may exist after 
the correction. As an example. Figure 7 contains two plots for Band 2. 
The top plot is a trace of the mean signal level of each of the first 136 
1 ijies in a radiometrically and geometrically corrected image (CCT-PT). 
The lower plot was constructed by computing the mean signal level for each 
detector (on a CCT-AT) over the full frame and using these mean values to 
construct eight scans (128 lines) of data. Note the similarity of the two 
plots and the 17-line periodicity of the top trace compared with the 



16-line period of the bottom. This 17- vs. 16-line period Is due to the 
resampling to a slightly smaller pixel size in the geometrically corrected 
data. Note that this periodic fluctuation in mean signal level Is not due 
to an error In the correction process; but stems from the round-off and/or 
truncation errors Introduced when applying calibration transformations to 
discrete signals. The result is low level striping effects which will 
appear In the final Image data even though the performance of the 
radiometric correction algorithm appears to be approaching its upper 
limits. 

Low Frequency Noise 

Analysis of low frequency noise (at or below scan mirror frequencies) 
produced the plots shown In Figures 8 and 9. In these figures, each trace 
Is a plot of the mean signal level as a function of scan number for a 
given detector. The mean signal level for a detector for a scan is 
defined as the mean of all the pixels In that scan line. Figure 8 
Illustrates the consistent difference between forward and reverse scan 
signal levels for Band 6. As suggested by Figure 3b, the mean signs! 
level for a forward scan is always lower than the mean signal for the 
following reverse scan, although the magnitude varies as one progresses 
down the frame. 

Band I, shown in Figure 9, contains several detectors which have 
significant low frequency variations. Detectors 4, 12, 10 and 8 have mean 
signal levels which vary between two states that are separated by 
approximately 2.0, 1.5, 1.0 and 0.75 signal counts, respectively. In 
other bands, this low frequency noise is much less evident, appearing in 
Band 2, Detector 1 ( 0.5 count), Band 3, Detectors 1 and 16 ( 0.5 count 
each), Band 5, Detector 10 ( 0.5 count) and Band 7, Detector 7 ( 0.5 
count). Note that this low frequency variation may remain in the higher 
or lower state for several seconds, and is not to be confused with the 
high frequency, within-scan noise reported previously (as for Band 7, 
Detector 7 )• [5]. This low frequency effect could be serious in 
water-related applications where average signal levels are low. 

Additional examination of the low frequency noise in Band 1 indicated 
that it was highly correlated between the noisy detectors (4, 12, 10 and 
8). One would expect that within a given scan, the signal variations of 
each detector would be best correlated with the variations of the adjacent 
detectors, and least correlated with detectors spaced further in the 
array. Although this is strictly the case for Band 4, and generally true 
with the other bands. Band 1 is markedly different. The inter-detector 
correlation matrices for all seven bands are presented in Figures 10-16. 
As well as the above mentioned correlation of detectors 4, 12, 10 and 8 in 
Band 1, the tendency towards every-other detector correlation (odd/even 
detector effect) is quite noticeable in Bands 2, 3 and 7. The cause of 
this low frequency noise is under investigation; a preliminary hypothesis 
is that it is the result of the DC restore procedure locking onto the high 
(low) phase of the coherent noise reported elsewhere [6]. 


SUMMARY AND CONCLUSIONS 


The overall quality of Thematic Mapper image data is very good. The 
spatial resolution is excellent and the radiometric correction procedures 
currently in use are working well, although some residual effects may 
remain which can cause low-level striping. Small but significant 
scan-direction effects exist in the reflective bands; they appear to be 
correctable in the two scenes examined by use of an empirical model 
employing exponential decay terms. Refinement and extension' of the 
initial correction model await further understanding of the sensor 
phenomena driving the effect. Band 6 has a unique scan-direction effect, 
markedly different from the effect observed in the reflective bands. 
Highly correlated low frequency noise in Band 1, Detectors 4, 8, 10 and 
12, may be due to a coherent-noise/DC-restore interaction; however, this 
phenomenon requires further investigation. It is recommended that any new 
radiometric correction procedures be integrated into a single process to 
minimize quantization and round-off effects. 

The shorter wavelength of Band 1 and a scan-angle range larger than 
previous Landsats both increase the magnitude of scan-angle effect due to 
atmospheric and scene bidirectional reflectance characteristics in 
Thematic Mapper data.' For a number of applications, approaches to reduce 
these effects should be developed. 
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CORRELATION BETWEEN DETECTORS OF LOW-FREQUENCY NOISE 

TM SAND 3 


Correlation Matrix for Scan Line MEANS <4> BAND: 4 CCT-BT 40049162G2 (IOWA SCENE) 
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FIGURE 13 

CORRELATION BETWEEN DETECTORS OF LOW-FREQUENCY NOISE 

TM BAND 4 


Correlation Matrix 'for Scan Line MEANS <5> BANDj5 CCT~BT 40049113262 (IOWA SCENE) 
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Correlation Matrix for Scan Line MEANS <7> BAND:7 CCT-BT 4004916262 {IOWA SCENE) 
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FIGURE 16 

CORRELATION BETWEEN DETECTORS OF LOW-FREQUENCY NOISE 

TM BAND 7 


